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ABSTRACT 

We present a V — I color-magnitude diagram for a region 1' — 2' from the center 
of M32 based on Hubble Space Telescope WFPC2 images. The broad color-luminosity 
distribution of red giants shows that the stellar population comprises stars with a wide 
range in metallicity. This distribution cannot be explained by a spread in age. The blue 
side of the giant branch rises to Mj ~ —4.0 and can be fitted with isochrones having 
[Fe/H] ~ —1.5. The red side consists of a heavily populated and dominant sequence 
that tops out at Mi ~ —3.2, and extends beyond V — I = 4. This sequence can be fitted 
with isochrones with —0.2 < [Fe/H] < +0.1, for ages running from 15 Gyr to 5 Gyr 
respectively. We do not find the optically bright asymptotic giant branch stars seen 
in previous ground-based work and argue that the majority of them were artifacts of 
crowding. Our results are consistent with the presence of the in/rarec?-luminous giants 
found in ground-based studies, though their existence cannot be directly confirmed by 
our data. The tip of the metal-poor portion of the giant branch occurs at the luminosity 
expected if M32 is at the same distance as M31 but is too sparsely sampled by this data 
set to provide a precise distance estimate. At fainter magnitudes, the rising giant 
branch is significantly wider (FWHMy_/ ~ 0.6 mag down to Mj ~ —1.0) than can be 
accounted for by photometric uncertainties, again due to a metallicity spread. There 
is little evidence for an extended or even a red horizontal branch, but we find a strong 
clump on the giant branch itself, as expected for the high metallicities inferred from the 
giant branch. If the age spread is not extreme, the distribution of metallicities in M32 
is considerably narrower than that of the closed-box model of chemical evolution, and 
also appears somewhat narrower than that of the solar neighborhood. Overall, the M32 
HST color-magnitude diagram is consistent with the average luminosity-weighted age 
of 8.5 Gyr and [Fe/H]~ —0.25 inferred from integrated spectral indices, extrapolated to 
the same radius and analyzed with the same population models. 

Subject headings: galaxies: abundances, galaxies: elliptical and lenticular, galaxies: 
evolution, galaxies: individual: NGC 221, galaxies: Local Group 
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1. Introduction. 

As the nearest example of an elliptical galaxy, M32 serves as an important beachhead in 
our campaign to understand the age and metallicity mixture of stellar populations in elliptical 
galaxies. While M32 may be somewhat unusual given its proximity to M31, spectroscopically it 
is identical to other faint ellipticals (Faber 1972; Burstein et al. 1984; Bica, AUoin, &: Schmidt 
1990; Gonzalez 1993; Trager et al. 1996). A key objective is to obtain deep color-magnitude (CM) 
diagrams to test conclusions thus far based solely on integrated colors and spectral index work. The 
proximity of M32, combined with the superb resolution afforded by the Hubble Space Telescope Wide 
Field/Planetary Camera 2 (WFPC2) permits a remarkable improvement in our ability to study its 
stellar content and brings us within range of this goal. 

In studying a stellar population, the basic variables we wish to determine are the mean stellar 
age and metallicity, and the distribution functions around these means^. For old stellar popula- 
tions like that in M32, it is difficult to measure both age and Z simultaneously from integrated 
light, as broad-band spectral shape and metal lines tend to vary together in lock step with both 
these variables (Worthey 1994). To break this "age-Z degeneracy", it is necessary to add another 
measured quantity that is specifically sensitive to age. Four such quantities have been suggested for 
M32: 1) Balmer line strength (Burstein et al. 1984, Faber et al. 1992), which is directly sensitive to 
turnoff temperature, 2) Sr IIA4077/Fe IA4045 (Rose 1985, 1994), which is sensitive to the ratio of 
dwarf to giant light and hence to turnoff luminosity, 3) ultraviolet surface-brightness fiuctuations 
at 2500A (Worthey 1993), again sensitive to turnoff luminosity, and 4) (1500-V) color (Bressan, 
Chiosi, & Tantalo 1996), which measures age through its effect on the number and luminosities of 
hot post-RGB stars. 

Many integrated light studies have claimed to have found evidence for a relatively young, 
few-Gyr-old stellar population in M32 {e.g. Baum 1959, Faber 1972, O'Connell 1980, Burstein 
et al. 1984, Rose 1985, Bica, Alloin & Schmidt 1990, Hardy et al. 1994). However, in retrospect it 
is clear that some of these investigations lacked the basic spectral information needed to distinguish 
age and Z, while others lacked sufficiently reliable population models to correctly interpret the 
important spectral features. 

Rapid progress is being made on both fronts. Stellar population models, taking both age and 
Z into account, have been constructed by Worthey (1994), Bressan et al. (1994), Buzzoni (1995), 
and by Tantalo et al. (1996). Accurate measurements of H/3 (Gonzales 1993, hereafter G93) and 
Sr/Fe (Rose 1994) in M32 are also now available (the former as a function of radius out to 1 Re). 
The upshot of this recent work is to confirm the earlier findings of relatively young stars, but the 
conclusion is now considerably more secure. 

Nevertheless, there are important assumptions in these population models that need to be 



^For some elliptical populations, the Fe-peak/a-element ratio also appears to vary (Worthey Faber, & Gonzales 
1992), but fortunately M32's element ratios seem close to solar, with Fe/Mg perhaps just slightly enhanced. 
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checked against a real CM diagram. For example, Worthey's models for the age and Z of M32 
assume a pure red clump for the horizontal branch (HB). If the H/3 were concentrated at intermedi- 
ate F-type temperatures, the need for a young population would be completely removed (Burstein 
et al. 1984). All models also assume a luminosity function for RGB and AGB stars that crucially 
affects the predicted red and infrared continua, spectral features, and surface-brightness fluctua- 
tions. Finally, the CM diagram permits measurement of the spreads about mean properties in a 
way that integrated light can never do. These spreads are as important as the mean values in 
decoding the star-formation history of the galaxy. 

This paper uses a new upper CM diagram of M32 from HST to address these and other funda- 
mental aspects of the stellar population. The major focus is on the metallicity and its distribution. 
That is because giant-branch loci are much more sensitive to Z than age, so conclusions about Z 
are relatively independent of the assumed age. The breadth of the giant branch implies a wide 
range in Z, from roughly solar down to below —1.0 dex. Paradoxically though, the distribution is 
rather narrow when judged by its FWHM compared to either the closed-box model of chemical 
enrichment or the distribution of metallicity in the solar neighborhood. Evidence for a wide giant 
branch and broad total range in Z was available before from pioneering ground-based CM diagrams 
in Fand / by Freedman (1989) and Davidge &: Jones (1992). However, it is clear now that the old 
data were just the tip of the iceberg — the present data give a much fuller picture of the upper 
giant branch, and a correspondingly greater total width 'm.V — I . 

A second focus is comparison with mean population age and Z deduced from ground-based 
integrated spectral indices. We can infer these values using Worthey models for the HST field by 
slight outward extrapolation of G93's data. The resulting age-Z pair, 8.5 Gyr and —0.25 dex, agrees 
perfectly with the color of the giant-branch stars in the CM diagram. This indicates a degree of 
consistency in the Worthey models but does not independently establish the age. Ages could be 
inferred from very deep observations of main sequence turnoff stars, but crowding precludes this in 
M32 inside about 3 arcmin, even with HST. Since upper giant-branch tracks are degenerate in age 
and Z much like integrated colors and metallic lines, it is not possible to derive an age from the 
upper CM diagram alone. Aside from verifying the red clump assumption of the Worthey models 
(in itself an important contribution), the new CM diagram has little to say about age. 

Nonetheless, the new observations represent an enormous gain over previous ground-based 
efforts in our ability to probe the stellar population of M32. They offer another tantalizing glimpse 
of how this so very interesting galaxy may have formed, and point the way for future work. 

2. Observations. 

Images of M32 were taken on October 22, 1994, with WFPC2. Two fields were imaged, one of 
M32 (POSl) and the other of a "background" field (P0S2) to determine the colors and magnitudes 
of the M31 stars that make up the majority of field stars at POSl. The POSl coordinates were 
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chosen so that the nucleus of M32 was roughly centered in WF3. The Planetary Camera (PC) 
image (0^'0455 per pixel) covers a region extending from 1 to 2 arcminutes almost due south of the 
center of M32 (see Figure | = Plate 1). The P0S2 field is situated ^ 5'.6 NE of POSl in an area 
originally intended to match the background surface brightness of M31 underlying M320. Four 500s 
exposures were taken in each of F555W ( V) and F814W (7) in each pointing. An additional three 
1000s exposures were taken through F300W, but we defer discussion of these data to a later paper. 
The coadded F814W image is shown in Figure § (Plate 2). 

Stellar photometry was carried out using the crowded- field, PSF-fitting package ALLFRAME 
(Stetson 1994), the most recent development in the DAOPHOT series of photometry packages. 
ALLFRAME differs from its predecessors primarily in its ability to use information from many 
individual frames simultaneously. The extra geometric and photometric information available from 
multiple frames extends the range of magnitudes and crowding conditions for which useful pho- 
tometry is obtainable. ALLFRAME is applied to the data after object detection and aperture 
photometry are carried out using DAOPHOT II (Stetson 1987). The global characteristics of the 
CM diagram produced using DAOPHOT II and ALLSTAR were essentially identical to those found 
using ALLFRAME. However, the latter produced fewer outliers and gave a cleaner-looking result. 

To push the star detection limit to as faint a level as possible, we coadded all images taken 
through the same filter. ALLFRAME was then applied to the F555W and F814W images simul- 
taneously. While analyzing the coadded images does not make full use of the information available 
in all 8 individual frames, experiments using both methods showed that our completeness fraction 
was significantly higher when coadding frames prior to object detection. 

Neither POSl nor P0S2 had a sufficient number of bright and reasonably isolated stars to allow 
proper characterization of the PSF. The PSF was consequently derived from F555W and F814W 
images of the globular cluster NGC 6397 observed with HST approximately two months prior to the 
date of the M32 observations. No focus changes were commanded during the intervening period, 
and desorption of volatiles will not have significantly altered the focus over this time. 

Owing to severe crowding, we ran three DAOPHOT II detection passes to identify and measure 
faint stars and to improve the photometry for the bright stars. After each pass, we subtracted all 
detected stars from the images and reran the finding algorithm to search for any newly-revealed 
faint stars. Three such passes yielded a candidate list of roughly 24,000 objects. Passing this list on 
to ALLFRAME resulted in the elimination of ~ 4000 objects, leaving a total of 19,969 measured 
stars. The final CM diagram is approximately one magnitude deeper at the 50% completeness level 
than it had been using a single detection pass. 

Identical procedures were used to find and measure stars in the P0S2 field, yielding a total 
of 4111 stars. Under ideal conditions, correcting for M31 background stars in the M32 field would 
require that, for each star found in P0S2, we subtract from the POSl sample the star most similar 



^Due to a pointing error, this match was not quite achieved — see below. 
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in color and magnitude. However, image crowding is clearly more important for POSl than P0S2, 
and proper accounting needs to be taken of their differing levels of completeness. Moreover, owing 
to a miscalculation, the P0S2 field as observed lies about 57" outwards from the M31 isophote 
that passes through POSl. The number of M31 stars we see in P0S2 consequently underestimates 
the number we would expect to find in POSl. Prom interpolation of blue M31 isophotes plotted by 
Hodge & Kennicutt (1982) and by deVaucouleurs (1958) and from color-index data by Walterbos 
& Kennicutt (1988), we estimate the F surface brightness gradient of M31 in the vicinity of M32 
to be approximately 0.0035 mag/arcsec. From this we conclude that the stellar surface density 
expected at POSl is 1.2 times that found at P0S2. 

Completeness tests were carried out by adding ^ 500 artificial stars, with colors and magnitudes 
selected from grid points in the CM diagram, to each of the POSl F555W and F814W images. 
The pixel locations of the added stars were randomly chosen and identical between F555W and 
F814W frames. The frames were then processed using DAOPHOT II and ALLFRAME in a manner 
identical to that applied to the original data. The results of these tests are shown in Figures ^ 
and ^, which show that the 50% completeness level for most POSl stars occurs at / ~ 24.9 and 
V ~ 25.4. Virtually identical tests were carried out for P0S2. 

The M32 CM diagram was corrected for contaminating M31 stars as follows. For a star at 
P0S2 of magnitude / and color V — I , a number of stars 

iV = 1.2 X Ci{I,V - I)/C2{I,V - I) 

of similar magnitude and color were subtracted from the POSl sample, where Ci and C2 denote 
completeness fractions computed at POSl and P0S2, respectively, and the factor 1.2 corrects for 
the misplacement of P0S2. In most cases is not an integer, and stars were subtracted once, 
twice, or not at all as determined by whether a randomly generated number fell above or below 
the fractional portion of N. In total, about 3600 stars were subtracted from the POSl sample to 
account for contamination by M31. 

The completeness tests also showed that the level of crowding leads to relatively large photo- 
metric uncertainties, ranging from ±0.19 mag RMS at / = 22.0 to ±0.46 mag RMS at / = 25.0. 
However, the uncertainties in Fand / are correlated, and the color uncertainties range from ±0.08 
mag RMS at / = 22.0 to ±0.32 mag RMS at / = 25.0. This is ihustrated in Figure |, where we have 
plotted the colors and magnitudes of artificial stars as returned by the DAOPHOT/ ALLFRAME 
processing sequence. 



3. Discussion 
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3.1. Morphology of the Color-Magnitude Diagram 

The CM diagrams for the POSl and P0S2 fields are shown in Figures |5| and ^. For stars bluer 
than V — I = 3 the magnitudes have been transformed from F555W and F814W to Johnson-Cousins 

V and / using the coefficients of Holtzman et al. (1995). The V — I color terms determined by 
Holtzman et al. are poorly constrained for very red stars. They are approximately zero at V — I = 
2.5 (the reddest color for which good measurements exist) and we have consequently adopted zero 
color terms for all stars with V—I > 3.0. The magnitudes have also been dereddened assuming E(i?- 
V) = 0.08 (Burstein & Heiles 1982) and using the absorptions tabulated by Holtzman et al. (1995) 
for stars of K5 spectral type. An electronic version of the photometry table is available on request 
from CJG. 

The most striking feature in Figure |5| is the obviously composite nature of the CM diagram — 
there is a wide spread in the colors of the giants. Particularly interesting is the strongly enhanced 
red sequence that bends over at / ?a 21.2, and extends to extremely red {V — I > A) colors. 
While Freedman (1989) and Davidge & Jones (1992) showed that the giant branch covers a wide 
range in color, this extended sequence was not visible in either previous work and represents a 
significant change in the picture of the M32 RGB. The present CM diagram is consistent, however, 
with Freedman's noted incompleteness at faint V magnitudes and her corresponding insensitivity 
to very red stars. This insensitivity is obvious when one considers the effect of a F-magnitude limit 
given the morphology apparent in the V vs V — I diagram, illustrated in Figure 0. In short, the 
new HST measurements reveal the full extent and morphology of the upper giant branch for the 
first time. Even with HST, however, there may be some incompleteness for very red stars with 

V — I > 3. For example, we find 14 starlike objects brighter than I = 22 that have no visible 
counterparts at all in the F555W image. Most of these are probably extremely red giants with 
V-I >4). 

Below about / = 22, the giant stars appear to merge into a rather fat giant branch of more or 
less uniform width down almost to the clump. Based on the completeness tests, the width of the 
giant branch above / ~ 24 is significantly larger than can be accounted for by photometric errors 
(the color uncertainty at / = 23 is 0.125 mag RMS, whereas the observed giant branch width at 
this point is 0.253 mag RMS — we return to this point below). The core-helium burning stars seem 
to be concentrated in a red clump. We see no indications of RR Lyraes or an extended reddish 
horizontal branch, even though the completeness in the region of the diagram where they should 
appear (/ 24.6, V — I ^ 0.5) is ~65%. The form of the fall-off in the density of stars below 
/ = 25 is due to a combination of completeness effects and photometric uncertainties. 

3.2. Crowding and the Brightest Giants 

We observe a significant number of stars to be brighter than the first ascent RGB tip. These 
may be AGB stars. However, their luminosities are significantly fainter than the AGB stars iden- 
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tified by Preedman (1989). Assuming a distance modulus (m — M)q = 24.43 (Ajhar et al. 1996), 
we expect the tip of the giant branch for older /more metal-poor populations to occur at / ~ 20.4. 
We find 23 stars with I < 20.4 in POSl and only one such star in P0S2, giving a surface density 
of AGB stars^ of 70 arcmin^^. Though the absence of such stars in P0S2 argues that these stars 
must be associated with M32, they do not appear to be as strongly concentrated towards the center 
of M32 as the integrated hght. Freedman finds 91 stars brighter than / = 20.4 in her inner field 
(a 1.42 arcmin^ region of her field which encompasses POSl — see Figure |l|), giving a very similar 
surface density of 64 arcmin"^. 

However, whereas our /-band luminosity function goes essentially to zero at / = 19.8, both 
Freedman (1989) and Davidge & Nieto (1992) see significant numbers of stars brighter than I = 19.8 
extending to / ~ 18. Since Freedman's inner field overlaps with our own, we can test directly the 
hypothesis that this inconsistency is due to differences in resolution. We therefore convolved our 
PC image with a Gaussian kernel of FWHM= 0''6 to match Freedman's CHFT resolution. We 
then integer-binned the data to a scale of 0"18 pix"^, reasonably close to the 0''2 pix~^ scale at 
the prime focus of the CFHT. Finally, we scaled the counts to match the exposure times used 
by Preedman, adding ~ 3 "background" counts per pixel as well as noise to match the gain and 
readout noise characteristics of the RCA chip on the CFHT. To characterize the degraded PSF 
in the standard way, the same process was applied to the NGC 6397 observations, constructing a 
PSF from 27 stars still visible in the resulting image. We applied two passes of the DAOPHOT II 
star-detection routine and used ALLSTAR to carry out the PSF-fitting photometry. The results 
yielded 147 stars detected in both the F555W and F814W frames in the region of overlap between 
the PC field of view and that of Freedman. Of these, 68 could be matched with the photometry 
of Freedman, yielding a mean 7-band magnitude difference of only 0.03 it 0.35 mag. The fact that 
we could match only half the stars in the field we attribute to differing completeness levels, which 
depend on rather small differences in the noise properties. 

The resulting CM diagram and luminosity function are remarkably similar to those obtained 
by Freedman. While the photometry of the original HST image yields only 3 stars brighter than 
/ = 20, smoothing the image results in a total of 19 stars measured to be brighter than this 
magnitude. Moreover, the mean color measured in the degraded image of the stars brighter than 
/ = 22 becomes bluer by 0.5 magnitudes, from V — I = 2.3 to V — I = 1.8. This strongly suggests 
that most of the bluer bright giants appearing in the ground-based CM diagrams of this region are 
in fact blends of fainter stars. 

The existence of the in/rared-luminous AGB stars found by Freedman (1992) and Elston &: 
Silva (1992) using J and ii' photometry appears to be compatible with our observations, although 
we cannot identify them uniquely with V and / alone. Freedman finds 125 stars with ii' < 18 in a 



We use the term AGB loosely here, meaning all stars brighter than the metal-poor RGB tip. Metal rich AGB 
stars may be fainter than / = 20.4, so the surface densities we derive will necessarily underestimate the true number 
of post-horizontal branch stars at all metallicities. 
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100" X 40" field at similar distance from the nucleus as our i/S'T observations; 85 of these JC-bright 
stars are also detected in J. Given our smaller image area, we would expect to have ~ 35 stars 
with < 18 in our PC observations, with half of these having K < 17. Most of these stars have 
J — K > 1.1, and would thus be expected to have I — K > 3.5 if the Milky Way bulge giants of 
Frogel & Whitford (1987) can be used as a guide. We would conclude that most of the i^-bright 
stars thus have I > 20, placing them comfortably below the bright tip of our /-band luminosity 
function. The Frogel & Whitford giants further show little correlation between V — I and I — K for 
I — K > 3.5, so the X-bright stars could fall anywhere within the broad V — I spread at the top of 
the /-band CM diagram. In contrast, the brightest and bluest stars in the Freedman sample might 
be expected to fall above the tip of the present /-band luminosity function; however, there are only 
8 stars with K < 17 and J — K < 1.1 in the Freedman sample, giving an expectation of just 2 
such stars for the PC field. The interesting problem posed by the Freedman (1992) and Elston 
& Silva (1992) data is that the majority of their /C-bright stars actually have J — K colors much 
redder than the Frogel & Whitford (1987) bulge sequence, making them presumably extremely red 
in / — /T, and perhaps also V — I. As noted above, there are 14 stars with / < 22 that have no 
V counterparts — perhaps these are stars that correspond to a significant fraction of the /T-bright 
sample, and for which it will be very interesting to obtain NICMOS observations. 



3.3. The Metallicity Distribution 

We argue that the composite nature of the giant branch in M32 is due to a wide range in 
metallicity. Gradients in the spectral index observations (see below) imply that the (light-weighted) 
mean age of the M32 population increases with radius, so there may well be a mixture of ages present 
in the PC field; however, a spread in age can account for only a modest amount of the spread in 
giant properties implied by the CM diagram morphology, as we now show. In Figure 0, we show 
a set of three isochrones (Worthey 1994) with ages 2, 5, and 15 Gyr. Since the three tracks can 
be made to overlap by making only modest metallicity changes among them, the age structure is 
impossible to deduce from the color distribution of the giant branch alone. The tradeoff between 
age and metallicity for any assumed RGB track is approximately d\og \t\ = — 1.7(ilog where 
the minus sign indicates that a younger age is counteracted by increased metallicity.[] On the other 
hand, since age effects are less important than Z, the color spread must primarily reflect the intrinsic 
metallicity distribution, with ambiguity in age causing only minor uncertainty. This is evident in 
Figure |^ where we show RGB tracks for clusters with diverse metallicities, and in Figure |8|, where 
we show theoretical RGB tracks covering a wide range in both metallicity and age. 

We can demonstrate the degree to which metallicity must predominate over age effects as 
follows. At Mi = -1.6, the width of the red giant branch is FWHM(y — / ) 0.6 mag, giving 



^Note that the behavior here is close the the behavior dlog|r| = —\.hd\og\Z\ found for integrated colors 
and metal-line strengths (Worthey 1994.) 
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a{V — I ) = 0.26 mag. The photometric uncertainty at this magnitude is a(V — I ) = 0.14 mag, 
yielding an intrinsic RGB width oia{V — I ) = 0.22 mag, or FWHM(1/ — / ) = 0.51 mag. Now from 
Figure ^ we see that FWHM(y — / ) for AZ = 1.2 dex is ~ 0.5 mag, giving a change in color with 
respect to metallicity of 0.42 mag/dex. From the scaling rule above, the age dependence therefore 
goes as 0.42 / 1.7 = 0.25 mag/dex. Assuming that the observed width must be the quadrature 
sum of both effects, we define 

At = FWHM in logr 

and 

Az = FWHM in logZ 

Thus we have 

[0.42Az]2 + [0.24A^]2 = 0.51^. 

Setting a maximum allowable age range of logr =1.0 dex (1.5 to 15 Gyr) and solving, we find a 
minimum allowable Az = 1.1 dex. The data clearly mandate a large range in metallicity irrespective 
of the range in stellar ages which may be present. 

The CM diagram implies that the distribution of metallicity is fairly smooth, but with many 
more metal-rich stars than metal-poor ones. This is illustrated in Figure |9|, where we compute the 
metallicity histogram by laying down a mesh of 10 Gyr isochrones on the CM diagram covering a 
wide range of [Fe/H], and then counting the stars between them brighter than Mj = —2.4, where 
the uncertainties are fairly small. A small correction for differences in RGB lifetime at different 
metallicities was also applied. No correction was applied for F-band incompleteness for the reddest 
giants, so the highest metallicity bins slightly underestimate the true number of giants by up to 
14 stars. One sees that the majority of stars lie along the red side of the RGB, with a decreasing 
tail toward blue, metal-poor RGB loci. There are very few stars (< 1%) more metal-poor than 
[Fe/H]~ —1.5. Although some of us initially had a visual impression of bimodality in the CM 
diagram, with a slight paucity of stars between the main "blue" and "red" RGBs, the [Fe/H] 
distribution shows no evidence of this. The illusion probably stems from the fact that RGB color 
begins to change rapidly with metallicity at about [Fe/H]= —0.5 and higher. The effect of assuming 
a population age of 15 Gyr is also shown; the histogram shifts only slightly to lower metallicity. 

In Figure ^ we compare the distribution of metal abundance found in M32 with closed-box 
models of chemical evolution. These simple, one-parameter models are one-zone, with no gas 
infall or outflow and zero metal content initially, and assume instantaneous recycling of heavy 
elements {e.g. Searle & Sargent 1972). The photometric uncertainties for the stars considered here 
are relatively small, and the models have not been smoothed. Quite striking is the fact that the 
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observed metallicity distribution is significantly narrower than the closed-box models would predict. 
This region of M32 appears at first glance to be more monometallic than even the solar cylinder. 
On the other hand, if the main ridge of stars in the CM diagram is comprised of populations with 
a spread in age, then the youthful stars may move to the next higher [Fe/H] bin, and it is possible 
that the real abundance distribution resembles that of the solar neighborhood more than it first 
appears in Figure |9|. For example, the metal abundance of the main red component is [Fe/H]= —0.2 
if one assumes an age of 15 Gyr, [Fe/H]= —0.07 for 8 Gyr, but [Fe/H]= 0.01 for 5 Gyr. Even so, 
the general shape of the distribution remains more or less constant with age, and the result that 
there are few low metallicity stars is secure. 

A similar metallicity distribution derived for the P0S2 M31 background field is also illustrated 
in Figure ^, renormalized by a factor of 5.7 so that the total number of stars is the same as that 
for the corrected POSl. Relative to M32, the M31 population has a peak at similar Z but a higher 
tail towards lower metallicities. The P0S2 M31 diagram resembles the solar neighborhood more 
closely, with the caveat that the age structure is completely unknown. The paucity of metal-poor 
stars relative to the closed-box model is also inferred for the nuclei of M31 and M32 (Worthey, 
Dorman, & Jones 1996), so it is probable that this condition exists at most radii in both M32 and 
M31. 

Worthey et al. (1996) argue that the simplest explanation for the lack of low-mass, metal- 
poor stars is that normal processes of chemical enrichment typically operate to produce fewer 
low-metallicity stars than the simplest closed-box model predicts. That is, modifications to the 
simplest model (like variable yield, spatially inhomogeneous enrichment, or variable IMF schemes 
that result in fewer metal poor stars) may be required to mimic chemical enrichment in the real 
universe. However, it is also possible that metal-poor stars could be present at larger radii. Large- 
radius storage of metal poor stars is a prediction of models of monolithic galaxy collapse (e.g. 
Larson 1975, Matteucci & Tornambe 1987, Arimoto &: Yoshii 1987) in which the metal poor stars 
are created at large radius and enrichment proceeds at smaller and smaller radii until most of the 
gas is consumed. The abundance distribution at a given radius predicted by these models is always 
narrower than the closed-box model (Larson 1975). 

The predicted number of metal-poor stars in the collapse model can be expressed as a fraction 
of total galaxy mass. In M32, the radius of the POSl field encloses about 75% of the total light, 
using the parameters of Kent (1987) and assuming a deVaucouleurs r^/^ profile. Conversion to 
mass requires a M/L ratio as a function of radius. Older ages and higher metallicities drive M/L 
up, so most galaxies are inferred to have a higher M/L in the nucleus. However, M32 is likely to 
be an exception to this general rule owing to the sizeable subpopulation of young stars resident in 



the nucleus (see Section 3^). Various population mixtures were modeled and compared to both 
the nucleus and the POSl positions, and the best-fitting M/Lb ratios were roughly M/Lb = 3 to 
4.5 on a scale where Galactic globulars have M/Lb = 2.7 for both M32 locations depending on 
the age-metallicity mixture. The amount of mass enclosed within the radius of POSl (which we 
have taken to be 75%) is consequently about 10% more uncertain than our estimate of the enclosed 
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light. 

In the closed-box model, about 10% of the stars in M32 should be more metal-poor than 
[Fe/H]= —1 assuming an approximately solar yield (the percentage would be higher if the yield is 
less than solar). Thus, if we stipulate a largely circular orbit distribution for metal-poor stars to 
account for the lack of such stars in our sample, there would still appear to be sufficient room to 
store them in the outer ~ 25% of the galaxy at R > 2' . HST observations are scheduled for a field 
in the outskirts of M32 that encloses about 92% of the light. If these observations uncover copious 
stars of [Fe/H]= —2.5 to —1.5, this would strongly support an outside-in dissipational origin for 
M32. If the observations find few such stars, then either they did not form in large numbers, as 
Worthey et al. suggest, or they were tidally stripped by close passages with M31, a possibility not 
without supporting evidence (Faber 1973; Nieto & Prugniel 1987). 

3.4. Clump Morphology 

Turning back to the CM diagram, the morphology of the core-helium burning stars in M32 
is a "clump" instead of a "blue," "extended," or "red" horizontal branch as commonly seen in 
Galactic (and M31) globular clusters with metallicity less than around —0.7. This confirms the 
conclusion of Rose (1994) that a clump, rather than red horizontal branch stars, must dominate 
the core- helium burning stars to produce the observed spectral line indices and broad-band colors. 
A "clump" morphology is expected for metal-rich populations (Hatzidimitriou 1991), so the fact 
that the M32 clump is strong is consistent with the giant branch color distribution. 

The extent to which the core-helium stars lie exclusively within a clump can in principle 
set upper limits on the numbers and ages of stars with metallicities between [Fe/H]= —0.8 and 
— 1.5. Superficially, the CM diagram is consistent with a pure-clump population. If a// core-helium 
burning stars are in the clump, then the metal-poor stars in M32 must be younger than ~10 Gyr, 
as populations this young or younger are neither observed (Stryker, Da Costa, & Mould 1985) 
nor predicted (e.g. Lee, Demarque, & Zinn 1994) to have RR Lyrae stars or stars intermediate in 
temperature between the instability strip and the clump. This picture is contrasted against one 
in which most of the stars in M32 are very old. In that case the more metal-poor components 
should display bluer horizontal branches like the "red" one of 47 Tuc or the "extended" one of M3 
(a theoretical approximation of which is shown in Figure ^) . 

Unfortunately, there are so few stars with [Fe/H]< —0.7, and the photometric errors are so 
large at the base of the CM diagram, that it is not possible to detect or rule out blue horizontal 
branch stars unambiguously. Using Figure ^ as a guide, and assuming a 15 Gyr population, there 
are 128 stars more metal-poor than [Fe/H]= —0.7 (out of 842 brighter than Mj = —2.4). The ratio 
of the number of bright giants corresponding to the Figure ^ cutoff to the number of horizontal 
branch stars in either M3 or 47 Tuc is about 0.5, so we would expect over 200 metal poor horizontal 
branch stars to be present under the all-15-Gyr hypothesis. To test whether we see these stars or 
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not, false stars were added to the observed CM diagram in either a 47-Tuc morphology or an M3 
morphology, with random errors approximately as observed. The conclusion from this exercise is 
that up to several hundred 47-Tuc-like HB stars could exist in the CM diagram, yet be hidden in 
the broad envelope of the more numerous true clump stars. Bluer M3-like HB stars could also exist, 
but probably less than 200 of them. In short, the present CM diagram cannot strongly constrain 
the upper age limits of the most metal-poor stars in M32. 



3.5. Age Information from Integrated Light 

Although the present CM diagram yields a clear picture of the metallicity distribution in M32, it 
does less well regarding age, especially of any younger components that might be present. Integrated 
colors provide some information on age if the metallicity distribution is approximately known, but 
the answer is intrinsically slippery because of model-to-model differences in the prediction of colors, 
which translate into a 35% scatter in age estimates (Chariot, Worthey, & Bressan 1996). The 
intrinsic power of the isochrones to delimit age is also not that much higher than integrated colors: 
AlogZ = — 1.7AlogT (see above) vs. AlogZ = — l.SAlogr (Worthey 1994). In addition, with 
integrated indices (including colors), we are presently limited to deriving an average age rather 
than a star formation history. 

Worthey (1994) has demonstrated that, to break the age-Z degeneracy, it is necessary to add 
information beyond colors and metal lines, for example, Balmer line strength. Using this approach, 
we will first attempt to derive an age under the simplifying assumption that the population is 
characterized by a single age. G93 took long-slit spectra to obtain line-strength indices on the Lick 
system, which can be directly compared to the Worthey (1994) model predictions. Over the POSl 
PC field of view, the G93 gradients (extrapolated from 45") were averaged, weighted by amount of 
light: 

Jarea/^rdA 

where < I > is the average index, fir is the r-band surface brightness as a function of position 
from Kent (1987), and A is the area over which the averaging was performed. Given the shallow 
color gradients in M32, r-band will adequately mimic V, where the indices are actually measured. 
The hght-weighted mean radius of the PC field is 1.8 Re, where G93 H/? = 1.92, Mg b =2.99, and 
< Fe >= 2.42. Worthey (1994) models predict that II/3 relative to the combined index [MgFe] = 
(Mg b X < Fe >)^/^ can give a simultaneous estimate for a mean age and a mean metallicity, with 
an age in this case of 8.5 Gyr and [Fe/H]= —0.25 for R = LSR^. 

This is noteworthy because if 8-Gyr isochrones are laid down on the POSl CM diagram and a 
metallicity distribution computed as in Figure^ the mean metallicity of the stars is [Fe/H]= —0.25. 
This demonstrates (at least) internal model coherence as shown by the close agreement in mean 
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metallicity estimates from star counts in the CM diagram on the one hand vs. Une strengths on 
the other. Within model errors, which can be substantial (Chariot, Worthey, &: Bressan 1996), the 
integrated colors of the M32 POSl field are also consistent with the same age and Z. 

Toward the nucleus of M32 the metallic indices of G93 get only slightly stronger, while H/3 
increases in strength dramatically (confirmed by Hardy ct al. 1994). Differentially within the 
models, the only way H/3 and metallic indices can both increase in strength is if the mean population 
towards the nucleus is simultaneously younger and more metal rich. The G93 indices at the nucleus 
indicate a mean age of ^4 Gyr and a metallicity just less than solar. To within the uncertainties, 
the predicted color at the nucleus is the same as that at 1.8 Re, despite a factor of two difference 
in predicted age because the age effect is balanced by an average abundance enhancement. 

Although age is the simplest explanation for the observed H/3 gradient and the fiat metallic 
feature strength and color gradients, two other alternatives should be considered: large numbers of 
either blue straggler stars or blue HB stars, if concentrated toward the center, could increase H/? 
strength as observed. However, both of these alternatives would make for bluer colors and weaker 
metallic features toward the center, which are not observed. Lacking a plausible mechanism to 
counter these effects, one would still require an age or metallicity gradient in M32. In the case 
of blue HB stars or young A-type main sequence stars, significant numbers of them appear to 
be ruled out in M32's nucleus by the Ca II index (Rose 1994; Worthey, Dorman, &; Jones 1996). 
Moreover, in contrast to M32, other elliptical galaxies generally show strongly increasing metal line 
strengths, reddening colors, and slightly weaker or constant H/3 toward the center. If blue stragglers 
or blue HB stars are responsible for M32's high central H/3, we must then explain why M32 behaves 
differently than other ellipticals. 

Despite these favorable checks, there remains little detailed information on the age structure 
in M32 because all the details are averaged into one or two indicators. Dropping the assumption 
of a single age for the POSl field, we constructed limited two-age models to try to quantify the 
uncertainty. The models begin with an assumed base age. Based on this age, a metallicity histogram 
is compiled from the CM diagram data. Then different age/metallicity components are substituted 
for the most populous bin in various proportions, and the results are compared with G93 indices. 
The constructed models are at all times consistent with the CM diagram for bright giants. 

The (merely illustrative, not definitive) result for a base age of 15 Gyr is that about 7% 
(by total mass) of a 2-Gyr population is needed to match the G93 indices. Interestingly, this is 
consistent with the results of a near-infrared survey by Silva k. Bothun (1996) of ellipticals with 
strong nuclear H/3 and blue central colors. For a base age of 10 Gyr, 5% of a 5-Gyr population 
would suffice. For a base age of 8 Gyr, the best match is obtained if there is 25% to 30% of an old, 
15 Gyr population present. From the way these populations balance, it is clear that over half of the 
stars at the POSl radius must be as old or older than 8 Gyr. All of these models match the data 
as well as the single-age 8.5 Gyr model. The conclusion is that there are many ways of balancing 
subcomponents so that the see-saw lands at a mean age of 8.5 Gyr. This sort of modeling, which 
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takes into account CM diagram information as well as integrated light indices, is in its infancy. 



3.6. The Tip of the Red Giant Branch 

The tip of the red giant branch (TRGB) in metal-poor populations can be used as a distance 
indicator. Old, metal-poor giant stars have peak absolute /-band magnitudes of Mj ~ —4 mag. Da 
Costa & Armandroff (1990) calibrated the magnitude of the TRGB for the range of metallicities 
(—2.2 <[Fe/H]< —0.7) defined by Galactic globular clusters. With this calibration, Lee, Freedman 
& Madore (1993) compared TRGB distances to 10 nearby galaxies to those from Cepheids and/or 
RR Lyrae stars; they found excellent correspondence at a level of ±10% RMS. Subsequent studies 
{e.g., Sakai, Madore, &: Freedman 1996a, b) confirm the excellent agreement between the Cepheid 
distances and the TRGB distances for metal-poor giants. A key property of the TRGB distance 
indicator is that M/ varies by less than ~ 0.1 mag for [Fe/H]< —0.7. 

Use of the TRGB distance indicator for M32 is problematic, given the high average metallicity 
of its population. To date, no empirical calibration for the TRGB exists for populations with 
metallicities higher than defined by the Galactic globular clusters. Based on the models of Worthey 
(1994), the expected behavior of the TRGB for more metal-rich systems is shown in Figure |lO[ Lee, 
Freedman & Madore (1993) conclude that the method could be applied to systems as long as the 
galaxies show an appreciable population of low-mass, resolved red giant branch stars with [Fe/H] 
< —0.7 dex. Figures ^, |l^, and |l^ underscore the importance of obtaining color information before 
blindly applying the TRGB method. The redder, metal-rich giants in M32 are almost a magnitude 
fainter than the bluer, metal-poor giants, and an uncritical application of this method to the metal- 
rich component could lead to a serious systematic error. The TRGB method works by passing an 
edge-detecting or "Sobel" filter over the /-band luminosity function that generates a strong signal 



at locations where the luminosity function has sharp discontinuities. As can be seen in Figure [li- 
the M32 luminosity function over all metallicities lacks a well-defined edge. Selecting just the metal 
poor stars (those with V — I < 1.8 mag) does produce a luminosity function with a sharper edge, 
but unfortunately in this case, there are too few such stars in the PC field to define the metal-poor 
TRGB accurately. Application of the edge detector to just the metal-poor giants yields a TRGB 
magnitude of m/ = 20.75, or M/ = —3.7 mag. This value has a large uncertainty, however, and 
its difference from the expected M/ = —4.0 (for the assumed distance modulus of 24.43) is not 
significant. 



3.7. Luminosity Function 



The /-band luminosity function (LF) we derive for M32 is shown in Figure 11, along with 
theoretical LFs for selected age/metallicity pairs. The model LFs were generated using a pure 
clump, as opposed to a red or blue horizontal branch. The mismatch between the data and the 
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model LFs in the clump region is likely due to a combination of (i) ill-determined completeness 
corrections at very faint magnitudes, (ii) the asymmetric nature of the photometric uncertainties 
in this region, which we have characterized as Gaussian for simplicity, (iii) uncertainties and a 
possible spread in the theoretically predicted clump luminosities, and {iv) the possibility that 
severe crowding is causing us to overestimate the luminosities of very faint stars in a manner 



analogous to that for AGB stars described in Section |3.2| . Of the four model LFs shown, the 
8 Gyr/[Fe/H]= —0.07 model is the best match to the observed discontinuity and bump in the 
region 20.5 < I < 22.0. However, the differences among the models are sufficiently small that a 
slight change in the assumed distance to M32 could easily negate this distinction. 

There are fairly substantial discontinuities in the LF at / ~ 20.5 and / ~ 21. The latter (which 
is the more significant) is clearly a consequence of the young/metal rich giant branch, which is 
essentially horizontal in Figure |5| at this magnitude. Assuming (m — M)o = 24.43, the discontinuity 
at I ^ 20.5 is very close to that expected for old/metal poor stars. Indeed, if we generate a LF using 
only metal-poor stars bluer than V — I = 1.8 (dashed histogram), the discontinuity is considerably 
enhanced. This is the tip of the giant branch used in the TRGB distance method. The cut again 
illustrates the importance of ensuring that the sample is not contaminated by stars more metal 
rich than [Fe/H]~ —0.7 when applying the TRGB method. As we noted in the previous section, 
the Sobel filter actually locates the edge of the metal-poor TRGB at a slightly fainter luminosity, 
but with little significance. 

To improve our sampling of the bright end of the luminosity function, we also analyzed the 
less-crowded half of chip WF4. Using a reduction procedure essentially identical to that used for 
the PC, this yielded another 20,000 stars (this field is slightly further from the center of M32). With 
allowances for the greater degree of undersampling and differing levels of completeness, the mor- 
phology of the CM diagram and the luminosity function obtained in the WF4 field are completely 
consistent with the results from the PC. 



4. Summary and Conclusions 

We have analyzed WFPC2 images of M32. Based on ALLFRAME crowded-field photometry 
we conclude the following: 

• The observed spread in color among giant branch stars requires a substantial range of metal- 
licity. Allowing for a large range in population age will have only a slight effect on the range 
of metallicity inferred. 

• The metallicity distribution is smooth and strongly skewed towards metal rich stars. The 
peak of the metallicity histogram occurs at —0.2 <[Fe/H]< +0.05 for assumed mean ages 
in the range 15 Gyr to 5 Gyr for the stellar population. A low metallicity tail extends to 
[Fe/H]Ri -1.5. 
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• The raw metallicity distribution appears somewhat narrower than that of the solar cylinder 

although, if there arc strong age admixtures, the distribution might resemble the local one 
more closely. In any case, the M32 distribution is much narrower than the closed-box model 
of chemical evolution. It remains an open question whether the missing metal-poor stars are 
absent entirely from M32 or whether they are stored at large radii. 

• The M31 abundance distribution in the outer disk sampled by P0S2 resembles the solar 
neighborhood if the stars there are mostly older than a few Gyr. 

• The width of the giant branch remains significant at magnitudes well below the tip, also 
supporting the conclusion that the population has a large dispersion in metallicity. 

• There is a strongly enhanced red clump at the bottom of the giant branch, consistent with 
Rose's (1994) analysis based on broad-band colors and spectral indices, as well as the metal- 
licity distribution implied by the giant branch morphology. At the same time, we cannot rule 
out a small minority population of extended-red or even blue horizontal branch stars, which 
would be expected to be present if the most metal-poor stars in M32 are older than 10 Gyr. 

• We do not see the very luminous blue AGB stars found in previous ground-based investi- 
gations. Experiments reveal that many of these stars are artifacts of image crowding. In 
contrast, a small population of the i^'-bright AGB candidates identified by Preedman (1992) 
and Elston & Silva (1992) may be present in our sample but cannot be identified by their 
corresponding / and V luminosities alone. 

• Whereas the metal poor giant branch extends to / —4.0, the metal rich population reaches 
only to /r^ —3.2. This underscores the importance of obtaining color information and selecting 
only stars with [Fe/H] < —0.7 before applying the TRGB method to the determination of 
distances. 

• The M32 CM diagram is consistent with the integrated line indices of Gonzalez (1993), 
extrapolated to the same radius. Interpreted with Worthey (1994) models, the indices give a 
mean age of 8.5 Gyr and a mean abundance of —0.25 dex. If 8-Gyr isochrones are overlaid 
on the CM diagram, the mean stellar abundance from counting stars as a function of color is 
also —0.25 dex. 

• The Gonzalez indices indicate a strong radial gradient in the mean stellar population with 
radius, the nucleus being much younger and somewhat more metal rich than the outer parts 
studied here. This is because the metallic indices are nearly fiat but slightly increasing toward 
the center and the broad-band colors nearly constant, whereas H/? increases strongly inwards. 

• Two-age models are also consistent with both the CM diagram and the Gonzalez (1993) 
indices, showing that there are many ways to mix ages and metallicities to match the ob- 
servations. Depending on the age of the bulk of the stars, modest subpopulations of quite 
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youthful stars are allowed. However at least half the stars must be at least as old as 8 Gyr. 
We are hopeful that future work will allow us to better constrain the options. 
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Fig. 1. — Digitized Sky Survey image of the field containing M32. The outlines indicate the WFPC2 
POSl and P0S2 fields and the Canada-Prance-Hawaii telescope CCD frame examined by Freedman 
(1989). The field shown in Figure |2| is indicated by the heavy white hne. The entire field shown 
here subtends 15' on a side. 

Fig. 2. — Linear stretch of the F814W PC frame. The lower-left corner of the field is 56 arcseconds 
from the center of M32. 

Fig. 3. — Completeness fraction as a function of / and F magnitudes sampled at 4 different colors. 
Filled symbols and solid lines correspond to the POSl field, while open symbols and dotted lines 
derive from completeness tests in P0S2. 

Fig. 4. — Colors and magnitudes of artificial stars returned in the course of completeness tests. 
The open circles indicate the grid of input colors and magnitudes, and the points show the colors 
and magnitudes returned after processing by DAOPHOT/ALLFRAME. 

Fig. 5. — ALLFRAME color-magnitude distribution of stars in the POSl (M32) field, corrected for 
background contamination using the color-magnitude distribution found in the P0S2 field. The 
magnitudes and colors shown have been transformed to V and / using a slightly altered transfor- 
mation (zero color terms for the very reddest stars) from that given by Holtzman et al. (1995). 

Fig. 6. — ALLFRAME color-magnitude distribution of stars in the P0S2 (background) field. 

Fig. 7. — The corrected color-magnitude diagram of stars in M32 as in Figure ^. Approximate 
V-band cutoffs are shown for these HST data (at V ~ 26.0) and for the Freedman (1989) data 
(at V ~ 23.2). Giant branch fiducials from Da Costa & Armandroff (1990) are shown for Galactic 
globular clusters M15 ([Fe/H]Ri -2.2), NGC 6752 ([Fe/H]Ri -1.5), and 47 Tuc ([Fe/H]Ri -0.8). 
Also shown is the giant branch of the old open cluster NGC 6791 ([Fe/Hjsi -1-0.2, age~8 Gyr 
from Garnavich et al. (1994)). The reddest point on this fiducial represents the reddest giant 
in the cluster, but it may not correspond to the helium flash point because the RGB is sparsely 
populated. Three Worthey (1994) isochrones are superimposed to illustrate the degeneracy of age 
and metallicity in RGB color. The solid line is (age,[Fe/H])= (15, —0.30), dotted (5, —0.06), dashed 
(2, -1-0.15). The predicted location of a metal-poor extended horizontal branch is shown. The line 
at Mi = —2.4 is the cut above which stars were counted for the histograms of Fig. |9|. The error 
bars are from artificial star tests and refer only to color V — I = 1.34. 

Fig. 8. — Worthey (1994) isochrones are shown on the POSl data for [Fe/H] between —1.2 and 0.0 
in 0.2-dex steps. More metal-rich isochrones are redder. Isochrones aged 8 Gyr are shown as dotted 
lines; isochrones aged 15 Gyr are shown as solid lines. The hypothetical locations of metal poor, 
extended horizontal branches are shown for both ages, the younger one nearly 0.2 mag brighter. 
Two examples of how bright AGB stars are predicted to lie in this diagram are shown extending 
beyond the 15-Gyr RGB isochrones of [Fe/H]= —0.6 and —0.4. The errors are for color V — I = 1.34 
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as in Figure 0. 

Fig. 9. — Counts of M32 stars caught between 10 Gyr isochrones of 1 different abundance. Stars 
brighter than Mi = —2.4 have small enough errors and the color change with abundance is large 
enough to make this a reliable method of estimating the abundance distribution of stars in M32. 
The raw counts between isochrones spaced at 0.2 dex intervals is shown as the coarse histogram 
with the heavy line. A small correction to compensate for the expected variable number of giants 
as a function of metallicity was incorporated. The light-line histogram assumes 15 Gyr isochrones 
rather than 10 Gyr isochrones. The dashed-line histogram shows the distribution for P0S2 and 10 
Gyr isochrones, and the counts have been multiplied by 5.7 to normalize the total number of counts 
to be same as for POSl. For comparison, two closed-box models of yield log[Zyieid/.^0] = —0.5 
and —0.2 are shown. The M32 distribution is much narrower than the closed-box models. Two 
derivations of the abundance distribution in the solar cylinder are also shown (Rana 1991, Pagel 
1989). The M32 histogram shown is even more narrow than the solar cylinder distribution, but 
if M32 contains many stars which are substantially younger than 10 Gyr, they will have higher 
abundance, and should be plotted (at most) one bin to the right. The true abundance distribution 
may therefore resemble that of the solar neighborhood quite closely, except for near [Fe/H]= —0.5, 
where M32 definitely has fewer stars than the solar neighborhood. The M32 sample of stars is 7 
times larger than the sample used for the solar neighborhood studies. 

Fig. 10. — Absolute /-magnitude of the tip of the red giant branch as a function of metallicity, from 
the models of Worthey (1994). 

Fig. 11. — Background-corrected /-band luminosity function of the POSl field before (solid his- 
togram) and after (filled circles) applying completeness corrections. The error bars reflect counting 
statistics only. The theoretical luminosity functions are shown after convolution with the estimated 
photometric uncertainties and application of the completeness function. The 15 Gyr models are 
shown as solid curves while the 8 Gyr models are shown as dotted curves. The models have been 
normalized using the number of stars found with 22.25 < / < 22.5. The dashed histogram shows 
the effect of counting only those stars with V — I < 1.8. 
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